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SOURCE CODE APPENDIX 



The algorithm is implemented within the runAsServer function of the DEM class 
and some other functions called within this function. These functions are provided 
5 below. 

The software code is presented in the Courier font. Code corresponding to the 
specific lines of the algorithm (see the section heading "Algorithm for Dynamic Slip 
Control") are presented in the Bold Courier font. Commentary is presented in 
10 the Italic Roman font. 



The DEM::runAsServer Function 

void 

15 DEM: : runAsServer ( ) 
{ 

Alert *alert; 
int selection; 
int i; 

20 

#ifndef EMBEDDED 
fd_set rfd; 
char error [256] ; 

25 DEM: : set Prompt ( ) ; 

DEM: :makeDataDirectories ( ) ; 
#endif 



30 The following high-lighted code fragment corresponds to line I of the algorithm 
(Le,x = x 0 ). 

// Performs user-specified application initializations 
applnit ( ) ; 

35 

/ / Insantiates and initializes user-defined components 
for server 

// initial state 
#ifndef EMBEDDED 
40 if (callAppInitComponents) 

#endif 

applnit Components ( ) ; 
#ifndef EMBEDDED 

45 

if (loadFile) 

DEM: : loadComponents (loadFile) ; 
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#endif // EMBEDDED 



5 The following high-lighted code fragment corresponds to line 2 of the algorithm 
[i.e> t = current_time() /. 

// Time variables are initialized with system time 
systemStartTime « currentApplicationTime = 
10 transitionRealTime = dem__time(); 



The following high-lighted code fragment corresponds to line 3 of the algorithm 
(Le., k= I). 

15 transitionld =1; 



The following high-lighted code fragment corresponds to line 4 of the algorithm 
(i.e., forever { 

20 

while (1) { 



The following high-lighted code fragment corresponds to line 5 of the algorithm 
25 [Le., A = next_event_tfme(A, x) ]. The function DEM::findProaction is described in 
a following section. 

// Finds the earliest proaction (if any) and stores 
it in 

30 // DEM: : scheduledTransition. It also assigns the 

right values to 

// DEM: sfiringDelay and DEM: : scheduledComponent 
DEM : : f indProac t ion ( ) ; 

35 

The following high-lighted code fragment corresponds to line 6 of the algorithm 

[i.e., x = set_interrupt_timer(A - (current_time() - 7 ))]. The 

DEM r.waitForlnterrupt functions are described in a following section. 

40 #ifdef WIN 

#ifndef EMBEDDED 

selection = DEM: : wait For Interrupt (&rf d) ; 
#else // EMBEDDED 

45 selection = DEM: :waitFor Interrupt () ; 

#endif 
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#else // WIN 

#ifndef EMBEDDED 
5 if (enterprise () ) 

selection = DEM: rwaitForlnterrupt (&rf d) ; 
else 

#endi£ // EMBEDDED 

selection = DEM: rwaitForlnterrupt ( ) ; 

#endif // WIN 



The following high-lighted code fragment corresponds to lines 7 and 9 of the 
algorithm [i.e., 7 next = current_time() and 7 = 7 next ]. The 
DEMr.updateTimeVariables function is described in a later section. 

// Updates time variables according to rt and sim 
options . 

DEM: :updateTime Variables (selection) ; 

The following high-lighted code fragment corresponds to line 8 of the algorithm 
[Le. 9 x = update^, ( 7 next - 7 )) ]. The update function is described in a later 
section. 

// If time has elapsed components are updated 
if (firingDelay > 0) 

for(i = 0; i < realtlmeComponents • size; i++) 
realtimeComponents . elements [i] - 
>update( firingDelay) ; 

The following high-lighted code fragment corresponds to line 10 of the algorithm 
[i.e. y x = compute(r, k y x) J. 

// The action associated with the scheduled 
transition is executed and 

// the return event is stored in event. 

Event * event = ( scheduledComponent- 
>*scheduledComponent-> 

a() [scheduledTransition] ) (alert) ; 

The following high-lighted code fragment corresponds to line 11 of the algorithm 
(Le., k = k + 1 ). 
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The following high-lighted code fragment corresponds to line 12 of the algorithm 
(i.e., } ). 
} 

> 



The DEM::findProaction Function 

void 

DEM : : f indProac t ion ( } 
{ 

#ifndef EMBEDDED 

char error [256 ] ; 
#endif 

firingDelay = -1; 

s chedul edComponen t = 0; 

scheduledTransition = -1; 

/* For the first component in a transient state */ 
while ( transientStateComponents . size > 0) { 

Component *c = transientStateComponents . elements [ 0 ] ; 

if (c->transientStates() [c->q] == 0) { 
transientStateComponents . remove (c) ; 
continue ; 

} 

int n = c->npro() / *p = c->pro(); 
double (Component: :**g) () = c->g(); 
Transition *t = c->t(); 
double d; 

/* ... for all proactions in current component starting in 
state q */ 

for (int j=0; j<n; j++) { 
if (c->q != t [p [ j ] ] . f rom) 
continue; 

d = <c->*g[p[j]]) (); 

if (d == 0) { 
scheduledComponent = c; 
scheduledTransition = p[j]; 
firingDelay = d; 

transientStateComponents . remove (c ) ; 
return; 
} 

} 
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/* current component is in transient state but 
no exiting proaction is enabled V 
#ifndef EMBEDDED 

sprintf (error, "DEM-f indProaction : Component %d/%ld is in 
transient state %d but no outgoing proaction is enabled. \n H , 
c->cid(), c->id, c->q) ; 

writeErrorLog (error) ; 
#endif 

transientStateComponents . remove (c) ; 
// Force c to go into its error state 
c->q = 0; 

} 

/* For all proactive components */ 

for(int i=0; i < proactiveComponents .size; i++) { 
Component *c = proactiveComponents . elements [ i ] ; 
int n = c->npro(), *p = c->pro(); 
double (Component :: **g) ( ) = c->g(); 
Transition *t = c->t(); 
double d; 

/* ... for all proactions in current component starting in 
state q */ 

for (int j=0; j<n; j++) { 
25 if (c->q != t [p[j] ] .from) 

continue; 

d = (c->*g[p[j] ] ) () ; 

// ... check if there is an enabled proaction that can 
happen before 
30 // firingDelay 

if (0 <= d ScSc (firingDelay < 0 || d < firingDelay)) { 
scheduledComponent = c; 
scheduledTransition = p[j]; 
firingDelay = d; 
35 } 

if (firingDelay == 0) 
break; 

} 

if (firingDelay ==0) 
40 break; 
} 

} 



10 



15 



20 



45 The DEM::waitForInterrupt Functions 

int 

DEM: : waitForlnterrupt ( ) 
{ 

int selection; 

50 

#ifndef WIN 
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struct timespec timeout; 
#endif 

#ifndef EMBEDDED 
5 char error [2 56]; 

#endif 

fd_set rfd; 
FD_ZERO(&rfd) ; 

if ( interruptAlerts . size > 0) 

//If the alerts from interrupt list is not empty- 
selection = 0; 
else if (firingDelay == 0) 

// If firingDelay is zero then scheduledTransition must be 
taken 

// immediately, 
selection = 0; 
else if (firingDelay > 0) { 

#ifndef EMBEDDED 

// if firingDelay is > 0 then scheduledTransition is 
scheduled in 

// the future: what to do depends on RT and SIM options, 
if (realtimeEnabled == 1) { 

if (simulationEnabled ==0) { 
#endif // EMBEDDED 

// If firingDelay > 0, rt on, sim off then wait with 
timeout 

// at (firingDelay - computationTime ) . 
double computationTime = dem_time() - 
currentApplicationTime; 

double timeToWait = 

(f iringDelay-computationTime) <0 ? 0 : ( f iringDelay- 
computationTime) ; 

#ifdef WIN 

Sleep ( (DWORD) ( timeToWait *1000 ) ) ; 

// Sleep doesn't return any values, so there's no way to 
tell if an 

// interrupt was received at this point. Since this 
function is used 

// only on embedded version of Teja for Windows we assume 
that no 

// interrupts were received, 
selection = 0; 

#else 

timeout . tv_sec = (long) timeToWait; 
timeout . tv_nsec = 

(long) ((timeToWait - timeout . tv_sec ) *1000000000 ) ; 
selection = (nanosleep (Sctimeout , 0) == -1) ? 1 : 0; 
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#endif // WIN 

#ifndef EMBEDDED 
} 

5 else { 

// If firingDelay > 0, rt on, sim on then 
scheduledTrans i t ion 

// is taken immediately. 
// 

10 selection = 0; 

} 

} 

else { 

// If firingDelay >0, rt off, (sim off or on) then no 
15 realtime 

// license is available on the system. 
scheduledTrans it ion 

// should happen in the future, but, since there's no rt 

// license, it's just ignored (an error message is 
20 printed) . 

writeErrorLog { "Realtime event scheduling option not in 
license . \ n H ) ; 

sprintf (error, 

"Proactive transition %d in %s %ld after %f seconds 

ignored. \n" , 

scheduledTrans it ion, 

ClassDescription . elements [ scheduledComponent - 

>cid()] 

->className, 
schedu 1 edComponen t - > id , 
firingDelay) ; 
writeErrorLog (error ) ; 

(void) select (FD_SETSIZE, &rfd, 0, 0, 0 ) ; 
selection = 1; 

} 

tendif // EMBEDDED 
} 

else { 

// If firingDelay < 0 then no proaction is enabled. Select 
without 

// timeout. 

(void) select (FD_SETSIZE, &rfd, 0, 0, De- 
selection = 1; 

} 

return selection; 

} 

int 

DEM: :waitFor Interrupt ( fd_set * rfd) 



25 



30 



35 



40 



45 
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{ 

int selection = 0; 
struct timeval timeout; 

#ifndef EMBEDDED 

char error [2 56]; v 

setRfd(rfd) ; 
#else 

FD_ZERO(rfd) ; 
#endif // EMBEDDED 

if (interruptAlerts . size > 0) 

// If the alerts from interrupt list is not empty 
selection = 0; 
else if (firingDelay == 0) 

// If the firingDelay is zero then scheduledTransi tion 
must be 

// taken immediately, 
selection = 0; 
else if (firingDelay > 0) { 

#ifndef EMBEDDED 

// if firingDelay is > 0 then scheduledTransi tion is 
scheduled in 

// the future: what to do depends on RT and SIM options, 
if (realtimeEnabled = = 1) { 

if (simulationEnabled == 0) { 
#endif // EMBEDDED 

// If firingDelay > 0, rt on, sim off then select with 
timeout 

//at (firingDelay - computationTime ) . 

double computationTime = dem_time ( ) - 
currentApplicationTime; 

double timeToWait = 

(f iringDelay-computationTime) <0 ? 0 : ( f iringDelay- 
computationTime) ; 

timeout . tv_sec = (long) timeToWait; 
timeout . tv__usec = 

(long) ((timeToWait - timeout . tv_sec )* 1000 0 00 ) ; 
selection = select (FD_SETSIZE, rfd, 0, 0, ^timeout); 

#ifndef EMBEDDED 
} 

else ( 

//If firingDelay > 0, rt on, sim on then 
scheduledTransi tion 

// is taken immediately, 
selection = 0; 
} 
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} 

else { 

// If firingDelay >0, rt off, (sim off or on) then no 
realtime 

5 // license is available on the system, 

s chedu 1 edTrans i t i on 

// should happen in the future, but, since there's no rt 
// license, it's just ignored (an error message is 
printed) . 

10 // Select without timeout. 

writeErrorLog ( "Realtime event scheduling option not in 
license . \n" ) ; 

sprintf (error, 

15 "Proactive transition %d in %s %ld after %f seconds 

ignored. \n" , 

scheduledTransition, 

ClassDescription . elements [scheduledComponent- 

>cid() ] 

20 ->className, 

s cheduledComponen t - > id , 
firingDelay) ; 
writeErrorLog ( error ) ; 

25 selection = select (FD_SETSIZE, rfd, 0, 0, 0); 

} 

#endif // EMBEDDED 
} 

30 else { 

// If firingDelay < 0 then no proaction is enabled. Select 
without 

// timeout. 

35 selection = select (FD_SETSIZE, rfd, 0, 0, 0); 

} 

return selection; 

40 } 



The DEM: mpdateTime Variables Function 

void 

45 DEM: : updateTimeVariables (int selection) 
{ 

#ifndef EMBEDDED 

if (simulationEnabled ==0) { 
50 #endif // EMBEDDED 

if (selection == 0) { 
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// Sim is off and timeout expired 
transitionRealTime = dem_time(); 

slip = transitionRealTime - (currentApplicationTime + 
f iringDelay) ; 

firingDelay = transitionRealTime - 
current Appl icat ionTime ; 

currentApplicationTime = transitionRealTime; 

} 

else { 

// Sim is off and alert or interrupt was received 
transitionRealTime = dem_time ( ) ; 
slip = 0; 

firingDelay = transitionRealTime - 
current Appl icat ionTime ; 

currentApplicationTime = transitionRealTime; 

} 

#ifndef EMBEDDED 
} 

else { 

if ((selection == 0) ( interruptAlerts • size ==0)) { 

// Sim is on, selection is 0 and no interrupt was 
received 

transitionRealTime = currentApplicationTime + 
firingDelay; 

slip = 0 ; 

currentApplicationTime = transitionRealTime; 

} 

else { 

// Sim is on and selection is <> 0 (an alert or an 
interrupt has 

// been received) was received 
transitionRealTime = currentApplicationTime; 
slip = 0; 
firingDelay = 0; 

} 

} 

#endif // EMBEDDED 
} 



The Component:: update Function 

void 

Component: : update (double elapsedTime) 
{ 

for(int i=0; i<this->ncs ( ) ; i++ ) 
x[i] += elapsedTime*xdot [i] ; 

} 
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